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CALCULATION  OF  THE  HYDRAULIC  FRICTION  OF  THE  COOLING  CHANNELS  OF 
GAS-TURBINES  BLADE 


L.  N.  Odivanov. 


Is  given  the  derivation  of  the  calculation  formulas  for 
determining  the  hydraulic  friction  of  the  cooling  channels  of  blades. 
In  formulas  is  considered  the  effect  of  the  preheating  of  air,  losses 
from  friction  and  centrifugal  backwater'  the  capacity  of  cooling 
channe Is. 

For  blades  with  radial  (longitudinal)  cooling  channels  are  given 
the  calculation  formulas,  suitable  for  a prog ramming^^^sm all  ETsVM  [ 
- digital  computer*].  For  blades  with  transverse  channels  is  given  an 
example  of  the  calculation  on  the  basis  of  the  selection  of  the  loss 
factors  in  handbook.  Calculation  data  satisfactorily  coincided  with 
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the  experimental  data  of  other  authors. 

In  gas  turbine  constructioij^was  received  vide  acceptance  of 


bl ad'^a ir- cooled l The  diagrams  of  some  of  them  are  given  in  Fig.  1 
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Daring  the  1-ayiptrt 


of  cooling  channels  it  is  necessary  correctly  to 


estimate  hydraulic  losses.  In  the  existing  aethcds  of  calculation 
[1],  etc.Ais  not  exaained  the  combined  effect  of  soae  special 
^_f**W4K=e^peculiarit ies  of  air  flea  in  the  cooling  channels  of  the  , 
blades:  the  effect  of  the  preheating  of  air  froa  the  heat^-^rpplied 
through  walls  froa  gas;  the  heat,  isolated  froa  friction;  the  effect 
of  centrifugal  forces,  etc.  This  leads  to  certain  arbitrariness  in 
the  calculations,  frequently -art  ia  pedes  the  account  of  the  specific 
conditions  of  the  flow  of  coolant  and  the  generalization  of  the  loss 
factors. 


In  work  [2]  is  experiaentally  investigated  the  effect  of  the 
factors  indicated  for  a rotor  with  the  tlades,  which  have  ■iUij.ii* 
deflector  [(Fig.  2)  . The  clearance  between  the  external  wall  and  the 
deflector  was  aain tain/withstood  because  of  point  punch-outs  on 
deflector. 


In  works  [3,  4 ] is  investigated  the 


of  blade  with 


thre  transverse  motion  of  air  coolant,  with  its  ejection  in  the  area 
cf  trailing  edge  to  the  side  of  trough.  Experimental  data  [4,  5], 
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etc.  confirm  the  possibility  of  applying  reference  data  [6]  with  some 
corrections. 

In  the  present  work  is  proposed  the  method  of  the  calculation  of 
hydraulic  friction  and  determination  of  the  flow  rate  of  cooling  air 
for  blades  with  the  longitudinal  and  cross  flow  of  coolant  with  the 
assigned/pj~e ncr iberl  geometry  of  channels  and  parameters  of  air 
coolant. 

Below  is  set  forth  the  derivation  of  the  calculation  formulas  on 
the  basis  of  the  equations  of  energy  and  continuity. 


L 


I 


i L 
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Fig.  2.  Suliy uUflic-  of  th*  cooled  blade  with  transverse  channels. 


Pa  ge  5 . 


1.  Equation  of  energy. 


The  eguation  of  energy  for  air  flow  in  cooling  channels  (in 

relative  motion)  is  r^eeeEA/ written  in  tbm  form: 


t Q 


II  “ MUX 


--  = 0.  (I) 


Here  Q - the  heat,  applied  to  air  fron  the  walls  of  the  heated 

to 


channel;  is  f,  relative  speed;  u - peripheral  speed;  T - the 
temperature  of  air  in  channels;  is  the  specific  heat  of  air. 


I*4-«res  and  "/Ji//"  designate  the  paraneters  of  air  coolant  at 
the  duct  inlet  and  at  egu^vtt/ yield  from  it  respectively  (Figs.  1,  2). 


According  to  the  eguation  of  the  1«t  law  of  the  theraody naaics 


dL,  + dQ  = c,,dT  -,p-  . 


(2) 


where  Zy  is  the  work  lost  in  friction;  p and  y - the  pressure  and 


X 
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air  density. 


In  the  process  of  the  expansion  of  air  its  preheating  through 
walls  from  gas,  from  friction  and  from  centrifugal  compression  let  us 


consider^  tit u\i gji'flie  single  average^index~^e-le^fttQ^boe^d  n,  i.e.,  let 
us  accept  /'/v"  -=  const.  , /cvjfca (S.)  ypcJU' 


Q-C'iT^-TJ  'r-n—RTn 


'-m) -]>"  « 


It  is  known  that 


bni  ^ 

I’  dL,  = L,  •=  rJjjiii 


(4) 


where  ®»«'"  is  -'Speed  without  taking  into  account  aC  hydraulic 

losses. 


Taking  into  account  expressions  (3)  and  (4)  equation  (1)  takes 
the  form  

J _ ( _/VuX  \ fl 


= 2 _ — or 
*"•"  n - \ K ' 


( l *U4X  \ 

K) 


' "Lx  - ul  4 «»•;.  <'»> 


Cndex  alectromesous  it  is  possible  to  determine  from  the  known 
relationsh4p/ra  tio 


r 
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Iv 

n — 1 T„ 


outlet  temperature  from  channel  from  the  equation 


T = T + ±T*  — — - — — 

' MUX  ' U~  1 k 

2 /? 

A — 1 
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For  determining  'ay.114X  froa  value  wm^.„  It  is  accessary  to  find  the 
factor  of  hydraulic  loss. 


2.  Factors  of  hydraulic  loss. 


For  blades  of  the  type  a,  b,  c (Fig.  1)  » duct  los sr'tfaa*:  are 
composed  of  the  following  forms:  input,  exit  and  loss  from  friction 
in  channel.  It  is  convenient  losses  to  real  outlet  velocity 

o 

(P* 

from  channel,  then  f-r  - »* — -p-  either  according  to  expression  (4)  we 
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will  obtain 


w~ 

*HX 


or 


w- 

MMX 


1 


fir, 


b-e 


) J 


*r  ti~ 

MWX 


+ *; 


1 + ':> 


(8) 


The  total  energy  losses  in  channel,  referred  to  the  unit  of 
lss,  will  be: 


-s  r^i 

J 2Dr 


cil  + E 


Hence 


di  + 


(9) 


where  ^ is  factor  of  local  lossf^  C - the  drag  coefficient  of  the 

sections  of  channel;  D,  - the  hydraulic  disaster  of  the  sections  of 
channel;  £ is  length  of  channel. 
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If  we  approxi aatel y accept  a linear  change  in  the  losses  from 
the  wall  friction  on  the  section  of  channel,  then  the  first  term  of 
expression  (9)  will  take  the  form: 


t = :»l  ( y i 

-Dr  „ \ «\1U,  / 2 n,  HUX 


Speed  change  on  section  from  output^ flld  to  the  inflection 


point  of  channed^negligibl^little,  i.  eT^jpcan 


W - TSL' 


Page  7. 


Taking  this  into  account  the  formula  for  determining  the  coefficient 
of  the  losses  of  the  channels,  depicted  on  Fig.  1,  will  take  the 
fora: 


t _ / tt’n  \V  ► •»  / \ 

'*  5w  + U— <•») 


The  velocity  ratio  is  determined  from  the  equation  of  the 
continuity 


/'  /•  T 

I nwx  1 MMX  * R 


tr_ 


(tl) 


From  equations  (6),  (7),  (8),  (10),  (11)  ic  determine^  and  7"iiwx. 
In  this  case  pressure  at  output/yicld- from  holes  can  be  accepted 
approximately  * that  which  is  being  chang  according  to  linear  law 
from  the  gas  pressure  on  periphery  ^in  at  .mMn tree  (into  lattice)  to  P-„ 
- at  output/yield. 


^ 


% 

A 
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^FOOTNOTE")1 . Hore  precise  according  to  data  [7]  /'„Ux  a little 

O' 

with  a decrease  in  the  clearance  6 and  with  an  increase  in 


the  dynamic  head  of  air  coolan tw^ENDFOGTNOTE. 


For  the  blades,  depicted  on  Fig.  2,  are  characteristic  still 
ether  forms  of  losses  whose  coefficients  are  enumerated  below: 


the  coefficient,  which  evaluates  loss  from  air  distribution  in 
central  channel 


coefficients,  the  evaluating  entry  loss'jinto  hole  0 (Fig.  2), 
losses ^s^sKock  of  flow  into  wall  (internal  surface  of  the  spout  of 
blade)  and  fron  the  rotation  of  flew  in  the  direction  of  transverse 
channels  i,  and  V’ 


the  coefficient,  which  considers  losses  from  the  mixing  of  flows 
upon  t h eTr^m  g an  d fron  sudden  e«pa  rsion  from  the 


channels  of  back  and  trough  (CM. 


Are  examined  below  recommendations 


the  selection  of  the 


a 
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loss  factors. 


On  the  periphery  of  blades  (Fig.  lb,  c)  the  air  flow  at 
output/yield  from  hole  depends  on  the  effect  of  the  flow  of  gas;  the 
mechanism  of  interaction  is  very  complex;  therefore  the  coefficient 
of  leaving  losses  :..u*  can  be  deterained  experimentally.  Tentative 
calculations  show  that  W'0-3  (the  apper  limit  corresponds  to  small 
radial  gap  6).  For  the  blades,  shown  in  fig.  2,  ’„uv  it  is  possible  to 
determine  by  Borda-Carnot 's  formula:  i,MI  = — ~-x  'j. 


Page  8. 


In  works  [3,  4]  the  total  losses  are  determined 


experimentally;  their  valueAmore  at  the  sealler  relative  clearances  h 

= h/d  between  the  leading  edge  of  the  deflector  and  the  internal 

surface  of  blade.  According  to  experiments  [3]  !„  t;»  = 4 + 5 with  h z 

- 

1,  while  according  to  the  data  [4]  0.7-1.  7 h s 2. 


The  satisfactory  agreement  of  these  data  is  obtained,  if 
are  selected  on  graphs  for  circular  exhaust  shaf t^minos  with  hood, 
upon  consideration  of  additional  entry  loss  into  hole  (account  of  the 


evenness  of  entrance)  [1,  6],  The  coefficients  of  losses 


at  the 


into  cooling  channels,  can  be  found  from  handbook  [6]. 
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The  values  of  drag  coefficients  lfc  transverse  channels, 

obtained  in  experiaents  [4],  satisfactorily  coincide  with  ^calculated 


[6].  The  drag  coefficients  of  radial  channels 


and  >„>JX  are 


deternined  from  formulas  for  ducts.  One  must  take  into  account  that 

the  friction  of  the  rotating  radial  charnels  is  more  than 



/un-rota  led/7 ones.  According  to  experiments  [5]  excess  e 


1 5o/o. 


These  coefficients  of  losses  and  recommendation  regarding  their 
selection  require  further  refinement  (especially  for  rotor  blades) 
with  the  accumulation  of  experimental  data. 


3.  Determination  of  the  preheating  of  coolant. 

cpdT*  =*  dQ  \ iiJu.  (12) 

The  equation  of  energy  in  the  complete  parameters,  the 
differential  form  and  in  relative  motion  has  the  form: 

^quantity  of  heat,  applied  to  coolant  (with  flow  rate  per 
second  <$£  G)  through  walls,  on  the  elementary  section  of  the  length 
cf  channel  dz  is  equal  to: 
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where 


dQ o “ OdQ  - k0  ( T'r  --  ! *)  di. 


- — - 


K,  flr 


p • * 


*“■  T; 


z - coordinate  in  the  direction  of  a radius;  6 is  wall  thickness;  a - 
heat-transfer  coefficient;  \ - the  theraal  conductivity  of  wall;  fj 
- the  periaeters  of  section/ee*-  of  wall. 


By  converting  equation  (12) , we  will  obtain 


(17 


* = r k0  + (<•>/)*  (r  4 z)  ' ±L  . 


J cp 
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At  the  average/aean  value  of  quantiti.es  T*t  the  definite 
integral  of  this  equation  Tra*ns®^fros  0 to  T takes  the  fora; 

i)jc  -n  ■■  f ■ «• 

dntg^nc^ 


where  r,  « , rn  - radius  of  the 


of  hole;  /*  - the 
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stagnation  temperature  of  coolant  at  the 

k = 


in  relative  motion; 


/OOTNOT^  l.  in  book  [8]  is  given  the  formula  for  determining  the 
preheating  of  coolant  in  nozzle  blades  (w  = 0) .CendFOOTNOtIL 


4.  Determination  of  the 


of  coolant. 


Outlet  velocity  from  hole  is 


expression  (8) 


converted 


w-  ■ 


~ „ 1 ^ / I!  1 »U\)  ' "(.MX  UK  "r  ‘'  ll  ^ I :M  _ ) 


,+-fe- 


The  drag  coefficient  of  channel  with  technical  roughness  with  Re 
numbers  > 2300  it  is  possible  to  determine  from  the  general-purpose 
formula  of  Colbrook  and  Hhite  [6]  taking  into  account  an  increase  by 
15o/o  in  the  data  [5]  the  rotation 
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The  calculation  in  channel  is  performed  according  to  sections 
<^&£-output/-y±erid  fro*  which  fro*  equations  (6),  (7),  (13),  (14)  and 

the  continuity 


by  the  eethod  successive  approximation  th^yare  detereined 

w T 

HMX  />  ‘ IIUX  /" 


*t  output/yi«4-d  froe  the  last/iot+er  section^^  is  calculated 
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P MMX  p Ptt 


’HMX  XNX i 


2iz i*C. 


k 

i*-l 


The  entire  calculation  is  fulfilled  also  by  the  aethod  of 
errroaeh/approxiaations  G before  agreeneai t^Pauxp  with  PHU* 
assigned/p  ce  scribed 


The  calculation  is  programmed  for 


coaputation  ETsV n 


A/a<V/ 


The  calculations  show  that  the  correctly  approximate  equality 


n -1 


* — 1 


(T.+  0,5Ar*) 


['-Or)] 


k 


C21) 


which  is  used  below. 

For  the  blades,  depicted  on  Fig.  2,  due  to  the  presence  of  many 
different  foras  of  friction  in  coolant  passage  the  calcu la tion^fcc? 
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A 


convenient  ly^conduc 
sections  along  channel. 


uti velnM"  totaling  hydraulic  losses  by 


If  we  write  equation  (5)  in  stagnation  parameters  and  -*§- 
consider  relationabip/rahio  (21),  then  we  will  obtain  expression  fcr 
determining  the  theoretical  total  pressure  P,1  cn  the  appropriate 
radius  of  the  blade 


k 


2j~r*<7'» r 0•5•s^  , 


(22) 


where  AT*  - the  preheating  of  air  on  section  from  the  entrance  into 
radial  channel  before  its  distribution  into  transverse  channels. 

Pa  ge  11. 


The  tenperature  of  stagnation  was  determined  from  converted 
izpression  (7) 

r*~r\  + \r*. 


The  actual  the  total  pressure  at  the 

•ction  will  be 


(23) 

each 
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I<n  handbook  [6]  the  coefficients  of  losses  € are  given  on  the 
basis  of  experiments  with  small  pressure  differentials  under 


isothermal  conditions.  In  long  channels  with 


preheating  of 


coolant  proceeds  a considerable  change  in  the  losses  and  aJt 
incidence/ drop  in  the  total  pressure  from  heat  supply.  Taking  this 


into  account  total  pitot  loss  (in  the  channels  of  const anfj'sect ion) 
approximately  » will  be 

V*  * (-J-  - l)-^  + + 1 

(for  channels  ' = 


"fOOTNOT^  *.  The  more  precise  value  of  the  total  pressure  at 
output/yieW  from  section  with  the  changing  cross  section  can  be 
found  from  equations  (14)  and  (20).  With  the  taken  u a p o-f  G 

is  determined  at  first  'W  and  then  /««■  IIDFOOTROTE^Z^ 


For  computation  2 it  is  convenient  to  a-ppiry  gas-dynamic 
functions.  Hith  tlfe'/a'ssigned/pgescribod  in  the  first  approximation. 


•eonsuspt ion  is  determined  the  function  <7(')  = 


and  from  tables  - 


~~:3r 
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-r-xRri 


The 

is  determined  Ap*. 


and  at  the  selected  value  6 


The  calculation  is  fulfilled  by  the  method ^successive 


dlsu 


sect  ion 


approximation^  changing  G before  agreement  in  exit 

/'.NX,  ^ 

Obtained  by  consecutive  subtraction^ Ap*  sections)  and  from 
the  expression 


ft's 

from  &- 


oVr~ 

y(>w)— -7- 


Page  12. 


In  this  case  pressure  at  output/y icl-d  from  hole  can  be  taken  as 
equal  to  the  pressure  of  gas  p2  at  output/y4e!4-  f rom  lattice  >. 


^FOOTNOTE*  1 . Work  [2]  shows  the  negligible  effect  of  ejecting  effect 
on  trailing  edge  (from  the  flow  of  gas)  on  the  G^noam pt  ion  of  air 
coolan t. /ENDFOOTND  T E> 
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The  proposed  calculation  nethod  Bakes  it  possible  to  determine 
coolant  distribution  according  to  span  cf  the  blade  (Fig.  2)  . 


Example  of  the  calculation. 


JoscJe*^-c 

For  Ofltma’the  method  of  calculation  and  selection  of  the  loss 


factors  was  determined  the  flow  rate  of  the  air  coolant  through  the 

channels  of  blades  [2]-  The  loss  factors  were  selected  from  handbook 

[6]: 

— 0,53  0,5f>  - for  the  holes  of  the  rim  of  disk  with  transition 

to  deflector; 

— 1 ,58  - for  a deflector  (taking  into  account  the 

separation  of  flow  on  holes  on  leading  edge)  ; 

-=4.76 ->-6,42  - for  the  slots  between  deflector  and  wall  of  blade 
(was  designed  by  formulas  for  lamellar-tubular  radiators)  ; 

ic  = U4; 

'/  m.x  “0,141  h-  0,1X)2; 

=0-77. 
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Fig.  3.  The  dependence  of  the 


ir  flow  rate  on  the 


relatioi 


ton  o 


pressures  and  temperature  of  the  blades: 
- according  to  experiments  [3]. 


- the  calculation;  O 


Page  13. 


The  preheating  of  air  AT*  from  the  sections  of  channel  was 
determined  from  the  equations  of  heat  exchange  and  energy  (13). 
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The  results  of  the  calculation  -°J£ 

I HUX  ' KMX 

and  €sh»«  the  temperature  of  blade  T , were  given  in  Pig.  3,  where  was 
visible  the  satisfactory  convergence  of  calculated  and  experimental 
data  [2  ]. 
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